Brownian motion model of random matrices revisited 
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We present a modified Brownian motion model for ran- 
dom matrices where the eigenvalues (or levels) of a random 
matrix evolve in "time" in such a way that they never cross 
each other's path. Also, owing to the exact integrability of 
the level dynamics, we incorporate long-time recurrences into 
the random walk problem underlying the Brownian motion. 
From this model, we derive the Coulomb interaction between 
the two eigenvalues. We further show that the Coulomb gas 
analogy fails if the confining potential, V(E) is a transcen- 
dental function such that there exist orthogonal polynomials 
with weighting function, exp[—/3V(E)], where f3 is a symme- 
try parameter. 
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Nearly all physical systems that are of significance are 
quite complicated due to their complex structure or in- 
teractions or a combination of the two. Having complete 
knowledge about these systems is unrealistic. As a re- 
sult, statistical analyses where we include as much salient 
information as is possible constitute (pen-)ultimate the- 
ories. The well-known examples are some of the most 
important branches of physics like statistical mechanics, 
random matrix theory, information theory, and so on jl] . 
Our interest here is in random matrix theory [Q where a 
complex or disordered quantum system is modelled in a 
realistic statistical way by drawing its Hamiltonian from 
an ensemble of random matrices. This mathematically 
elegant theory becomes physical also as we incorporate 
dynamical features. The basic theme is to develop a 
probabilistic (statistical) framework which is systemat- 
ically augmented by Newtonian or dynamical (determin- 
istic) characteristics that are fundamental attributes of 
the system. In this Letter, we present a modified Brow- 
nian motion model, originally proposed by Dyson ||,Q, 
that takes into account basic features of level dynamics. 

Matrix models are now employed to understand phys- 
ical phenomena in varied fields, particularly for heavy 
nuclei [||, mesoscopic disordered conductors [||, two- 
dimensional quantum gravity Q, and chaotic quantum 
systems (8| . A lot of excitement and activity in this field 
has been going on for nearly two decades due to nearly- 
perfect modelling of universal conductance fluctuations 
in metals Q and local spectral fluctuations in chaotic 
quantum systems ||. 

Given an N x N random matrix with eigenvalues, 
Ei, E 2 , En, the joint probabilty distribution function 



P(E u E 2 ,...,E N )=C N0 e^ w ^ E ^-' E ^, (la) 

N 

W(E U ...,E N ) ='£v(E j ) -^loglA - Ej\ (lb) 



gives the probability for eigenvalues to be lying in unit 
intervals around the points, EiS. This is also the proba- 
bility density of N unit charges moving freely on the real 
line under the influence of forces derived from the poten- 
tial energy, (lb). Dyson HH generalised the notion of 
matrix ensemble so that (la) is the equilibrium distribu- 
tion of a statistical mechanical system. Here the many- 
body system has fictitious particles as levels or eigenval- 
ues and it evolves in some fictitious time. In general, the 
dynamical system is in some non-equilibrium state. This 
gives birth to the Brownian motion model where eigen- 
values, Ei perform Brownian motion with no well-defined 
velocities or inertia. The Coulomb gas is thus not a con- 
servative system, the particles are actually experiencing 
frictional forces that resist their motion. 

Two decades later, in an attempt to provide dynamical 
origin of (la), it was shown by Pechukas (!(]] and Yukawa 
flli"| that the levels can be thought to evolve with a pa- 
rameter, t (like "time") if one considers the Hamiltonian 
H = Hq + tV where V is Hermitian. It was shown that 
for eigenvalues, equations of motion can be written and 
the Hamiltonian underlying level dynamics is, in fact, 
an exactly integrable generalised Calogero-Moser system. 
There are two basic features that emerge out of these re- 
cent works : (1) due to the fact that levels evolve adia- 
batically, there should be no curves crossing each other 
as the system evolves and the many-body system equi- 
librates; (2) since the system is an A-body Hamiltonian 
system, in principle, Poincare recurrences are bound to 
occur. Regarding the recurrences, it is quite clear that as 
N becomes large, the recurrence time will be very large 
too. It is also quite obvious that these two important fea- 
tures are not satisfied in the Brownian motion model of 
Dyson. The aim of this Letter is to present a Brownian 
motion model where these dynamical inputs are incorpo- 
rated, in consonance with the opening discussion. 

Consider N levels or eigenvalues as random walkers 
which are initially close together, walking on a line with- 
out crossing each other, with a non-zero probability of re- 
currence after a long time. In this random walk problem, 
we wish to find the probability distribution of positions 
of the random walkers at some later time, t, given that 
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there is recurrence at time, 2t. At time t, the position 
of an eigenvalue is Ei(t), i — 1,2, N. Initially, let the 
positions be £";(()) = E{ t o satisfying 



— Emax < Eifi < £3,0 < ••■ < En,o < £max 



(2) 



where -Emax is a fixed bound. The objective is to find the 
joint probability distribution function for levels evolving 
without crossing each other's path, i.e., maintaining the 
inequalities : 

Ei-i(t') < Ei(t'), i = 2,...,N, for t'<t. (3) 

Let us denote the set of initial (final) set of levels, 
(E lfi ,...,E N<0 ) ((E 1} ...,E N )) by E Q (E). Similar ran- 
dom walk problem is encountered in the discussion on 
commensurate melting where domain walls move [ p"2"| . 
We use the arguments presented in this work below to 
our purpose. One can systematically treat this random 
walk problem on a lattice with discrete time. Since we 
are interested in long times, we take recourse to the con- 
tinuum approximation where we have the simple case of 
Brownian motion. The probability distribution of N free, 
unrestricted random walkers is |13| 



P° N (E ,E;t) 



(2irDt) N / 2 



(4) 



where in consistency with the Brownian motion, D is the 
diffusion constant that sets the time scale via 



((Ei - E h0 ) 2 ) = Dt. 



(5) 



The angular brackets denote the averaging with proba- 
bility distribution. The probability distribution satisfies 
the TV-dimensional diffusion equation, 



N 

E 



d 2 p 



N 



2 dP, 



N 



dx 2 



D dt 



(0) 



The distribution (4) must evolve in such a way that the 
linear manifolds, 



Ei — E, 



i+l) 



»=1,2,...,JV-1, 



(7) 



are never crossed. Equivalently, the distribution Pm = 
on these manifolds. Since the manifolds (7) are linear, 
the problem can be solved by the method of images. The 
levels are allowed to walk unrestrictedly and randomly 
with positive and negative weights starting at Eq and its 
mirror images, according to (4). If the initial weights are 
chosen such that Ppj is antisymmetric under reflection in 
each manifold in (7), all subsequent distributions will be 
antisymmetric too. 

To obtain an analytic expression, we use the elements, 
g of iV-dimensional symmetric group, Sjy Ej ], to permute 
the components of any A-tuple, E = {Ei, E 2l En}- 



To obtain Pjv, as remarked above, we start a level 
with weight, e(g) at each mirror-image point gEg, where 
e (.9) = +1 (~ 1) f° r even (odd) permutations correspond- 
ing to even (odd) number of reflections. For E satisfying 
(3), we have 

P N (E ,E;t)= J2 e(g)P N (9Eo,E;t) (8a) 

-|E| 2 -|gp| 2 
p (2-Dt) 

(8b) 



= A 



N (2irDt) N / 2 ' 
where ,4^ is the antisymmctriscd sum, 

A N (E Q ,E;t) = <9)e^ 
g&S N 



(9) 



Let us now define the recurrence of fictitious particles or 
the reunion of the walkers around E by 



E — E\ 



< Ei < E 2 < ... < E N < E + E r 



(10) 



where we assume (without loss of generality) that the 
initial and final mean positions are the origin and E. 
Thus, E.gE = £f =1 (Ej - E)(gE )j - O(^nrax)- Due 
to antisymmctrisation, the form taken by A n contains 
the Vandermonde determinants in initial and final com- 
ponents. That is, 



TV 



AslEr,E:l)= J*)~ nN J [\E : 



1!2L.(JV-1)! 



.,-) I ^max 
" Dt 



3,0 



Ek,o\\Ej 



j<k 



Ei, 



(11) 



n N — 7}N(N — 1). The asymptotic probability for recur- 
rence is, in fact, from (8b) and (11), given by 



P N (E ,E;t) 



e 2nt 



N 



{2v) N / 2 {Dt) N2 /- 

1 + °( : % X ) 
1!2!...(AT - 1)! ' 



(12) 



Writing (12) as the partition function for N fictitious par- 
ticles, we obtain the Coulomb interaction. To emphasise 
this important result, we have shown how the Coulomb 
gas analogy follows directly from a conservative (neigh- 
bors being always the same), recurring Brownian motion 
model. 

Any two eigenvalues interact with each other by the 
Coulomb interaction in two dimensions. The potential 
that confines these eigenvalues, V(E) is the one we re- 
quire now to complete our discussion. The universality 
of local spectral fluctuations is due to the fact that the 
Coulomb gas analogy holds true. The non-universal as- 
pect is the average level density, cr(E), which is connected 
with the confining potential by the following equation M : 



2 



dE 



, cr(E') _ dV{E) 
E-E' ~ dE ' 



(13) 



The left hand side is a principal value integral. 

With the confining potential thus found, one can ex- 
plore the situations where for the weighting function, 
exp[— (3V(E)\, there exist orthogonal polynomials. This 
enables one to use the powerful method of orthogonal 
polynomials in finding various correlations between the 
eigenvalues Q . We now come to our second result in the 
form of a theorem. 

Theorem : If the matrix elements of a random matrix are 
distributed by a transcendental function so that the con- 
fining potential of eigenvalues in the continuum approxi- 
mation is also a transcendental function V(E), then the 
Coulomb gas model breaks down if there exist orthogonal 
polynomials with a weighting function exp[— (3V(E)\. 
Proof : Using Reductio ad absurdum, let us assume 
that the Coulomb gas model works. Let the points, 
Ei, E2, -B3, ... be such as to minimize the potential en- 
ergy, 



A' 



i<j 



Then, 



dV 
dE { 



E,; — E.; 



(14) 



(15) 



Let us assume that there exist orthogonal polynomials, 
Pn{E), with e~P y ( E ^ as the weighting function. Let Ej 
be the zeros of p^ (E) . By the fundamental theorem of 
algebra, 

PN (E) = C{E - Ei)(E - E 2 )(E - E 3 ) (E - E N ), 

(16) 

C being a constant. Next, we can write Q 



E 



Ei — Ei 



(17) 



Substituting in (15), we get a differential equation for 
Pn(E), 



dV 



2p' N (E 3 ) 8E 



(Ej) = 0, 



(18) 



which is absurd for a ratio of two polynomials can never 
be a transcendental function. This concludes the proof. 

The Brownian motion model has served as a basis for 
most of the developments in quantum transport theory. 
In this Letter, we have pointed out the features that are 
now well-known from the results in level dynamics. Since 
Dyson's Brownian motion model does not satisfy these 
requirements, we have given a modified model. This has 



the most important distinction of giving us the Coulomb 
gas model as a result of the constraint that the levels 
cannot cross each other. Thus the Fermionic nature of 
levels emerges from the random walk problem considered. 

Moreover, we have proved that the Coulomb gas anal- 
ogy breaks down for the case when the confining poten- 
tial is a transcendental function such that there exist or- 
thogonal polynomials with the corresponding weighting 
function. However, even if the Coulomb gas model is 
inconsistent with any transcendental confinement, if the 
joint probability distribution function (j.p.d.f.) satisfies 
the iV-dimensional diffusion equation, the local correla- 
tions may still be universal. This situation arises in meso- 
scopic disordered conductors Jl^] where the conductance 
fluctuations are universal and the confining potential is 
log 2 (_E). The Coulomb gas analogy breaks down there, 
however the j.p.d.f. satisfies the diffusion equation. 
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